34 International Conference on High Energy Physics, Philadelphia, 2008 

Mixing and CP Violation in Decays of Charm Mesons 

B. Meadows * 
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The phenomenon of mixing in neutral meson systems has now been observed in all flavours, but only in the past 
year in the D° system. The standard model anticipated that, for the charm sector, the mixing rate would be small, 
and also that CP violation, either in mixing or in direct decay, would be below the present levels of observability. It 
is hoped that further study of these phenomena might reveal signs of new physics. A review of recently available, 
experimental results is given. 

1. INTRODUCTION 

Mixing and CP violation (CPV) in the neutral D system were first discussed over thirty years ago [1] but mixing 
was observed for the first time only very recently [2, 3]. Since then, these observations have been confirmed [4-7] in 
other experiments and in other D° decay modes. 

Unlike the K, B and B s systems, for which mixing was observed years earlier, the short distance (AC = 2) 
amplitude contributing to mixing in the D system arises from box diagrams with down- rather than up-type quarks 
in the loops. The d and s components are GIM-suppressed, and the b component is suppressed by the small V u b 
CKM coupling. In the standard model (SM), therefore, long range, non-perturbative effects, a coherent sum over 
intermediate states accessible to both D° and D°, are the main contribution to mixing. These are hard to compute 
reliably, however [8-11]. 

1.1. Notation and Formalism 

D° and D° mesons are produced in flavour eigenstates, but propagate in time t and decay as mixtures of eigenstates 
D\ and Z?2, with masses and widths M\fl and Ti,2, related to the flavour states by: 

|L>1> = p\D°>+q\W> ; \D 2 > = p\D°> -q\W> (1) 

It is usual to define mixing parameters x, y, r M and 4> M and a decay parameter A/: 

in which M and T are averages for the two mass eigenstates and Sf, <pf are strong and weak relative phases, 
respectively, for the amplitudes Af(Af) for the decays of D°(D°) to the final state /. In the absence of CPV in 
mixing p = q = 1/^/2 and 4>m (the mixing phase) is zero, D\ is CP-even and Di is CP-odd. If there is no direct 
CPV in decay, then either <f>f or Sf = Q. 

Recent estimates for |x| ~ 1%, \y\ ~ 1% are in good agreement with the present observations and are orders of 
magnitude smaller than values measured for the other neutral mesons. Since charm decays are dominated by the 
2x2 sector of the CKM matrix, where weak phases are absent, CPV, in either mixing, or in time-integrated decay 
rates, is unexpected. Either a large value for |jc| (e.g. > \y\) or the observation, at current sensitivities, of CPV in 
charm decays would be regarded as an indication of new physics [9-12]. 

In the following section, we discuss the experimental evidence for mixing as observed in time-dependent effects in 
hadronic decays. Then, in section 2.2 we present results from its effect on lifetimes. In section 4, we outline the 
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present state of experimental searches forCPV in time-integrated decay rates of charmed mesons in section 4. In 
section 5, we present a summary of the current knowledge on mixing parameters and outline possibilities for future 
progress. 



2. EVIDENCE FOR MIXING IN TIME-DEPENDENCE OF DECAYS 



It follows from Eqs. (1) and (2) that mixing gives rise to a time-dependence in amplitudes A{A)f^ for decays 
D°(D°) — > /(/). When the final state / is accessible to decay by both D° and D°, interference between direct decay 
and decay preceded by mixing will occur. Neglecting CPV and terms in x or y beyond second order, the decay rate 
for D — > /, for example, deviates from exponential 
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The first term in square parentheses is the direct decay rate, the third is proportional to the mixing rate {x 2 +y 2 )/2 
and the middle term comes from the interference. The latter makes it possible to measure the small quantities x and 
y to first order, but only if the strong phase difference Sf is known. 



2.1. Evidence for Mixing In D° — ► K + n Decays 

Observation of the deviations from exponential in the "wrong-sign" (WS) decays D° — > K + ir~ l , a signal for 
D° — D° mixing, has been attempted many times [13-18] without success. These very rare decays are Doubly- 
Cabibbo-suppressed (DCS) and must be compared, experimentally, with the far more copious, "right-sign" (RS) 
Cabibbo- Favoured (CF) decays to K~ix + . in which the deviations are negligibly small. 

The BABAR collaboration, using a huge 384/6~ 1 e + e~ sample of D° mesons, whose flavour was determined from 
the sign of the pion from the D* + — > D°tt + decays from which they came, were finally able to observe mixing. In 
Fig. 1(a), the ratio of WS to RS decays of the combined sample of D° and D° decays is shown as a function of 
proper time t. The average value of this ratio, ~ 0.36%, is clearly not constant. The mixing fit, including a linear 
rise corresponding to the second term in Eq. (3), is preferred. 




t igure 1: Evidence for D° - D° mixing. The ratio of WS /RS yields as a function of proper time t is shown for the (a) BABAR 
and (b) CDF experiments. The dashed plot in (a) shows the projection of the mixing fit onto the time axis, and has a \ 2 of 
1.5. The flat line (no mixing fit) has x 2 of 24. The dashed line in (b) is a parabolic fit giving values for the DCS rate, x' 2 
and y . In (c), the Heavy Flavour Averaging Group (HFAG) [19] summary of measurements of yep is presented. The mean, 
1.072 ± 0.257%, differs significantly from the no-mixing expectation of zero. 



1 Hereafter, unless otherwise stated, it is implied that charge conjugate modes are included. 
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The strong phase Sk* was poorly known, so this measurement was only able to determine values for y' = (9.7 ± 
4.4 ± 3.1) x 10~ 3 , x' 2 = (-0.22 ± 0.30 ± 0.21) x 10~ 3 and R D = (0.303 ± 0.019)%, where 



and R D is the ratio of DCS ro CF decay rates. The BABAR result was found to be inconsistent with a no mixing 
hypothesis, with a significance estimated at 3.9cr. These results assumed no CPV. Separate fits to D° and to D° 
showed no significant difference. It was concluded that no CPV was observable at the current precision. 

Using a similar technique, but with data from 1.96 TeV pp collisions, the CDF experiment has confirmed this result 
with a 3.8cr significance [5]. The plot of WS/RS ratio vs. t is shown in Fig. 1(b). Assuming no CPV, they obtained 
values x' 2 = (-0.12 ± 0.35) x 10~ 3 , y' = (8.5 ± 7.6) x 10~ 3 and R D = (0.304 ± 0.055)%, in excellent agreement with 
those from BABAR. 



Deviations from exponential decay, introduced by mixing, alters the effective lifetimes t for D° decays in ways 
that depend upon the CP symmetry of the final state [20, 21]. Decays to K~n + (mixed CP) and to K~K + or ir~ii + 
(CP-even) have been used in searches for mixing by measuring the quantities 



Both are zero in the absence of mixing or CPV . Clear evidence for mixing (y C p ^ 0) has now been seen in this way 
by both the Belle [3] and BABAR [4] collaborations with significances of 3.2cr and 3. Oct, respectively These results, 
with those of earlier, unsuccessful searches for mixing, are summarized in Fig. 1. 

No evidence for CPV (A T ^ 0) is yet observed, however. Belle measured A T = (0.010 ±0.300 ±0.150)% and BABAR 
obtained 2 A T = (0.260 ± 0.360 ± 0.080)% giving an average A T = (0.123 ± 0.248)%. 

2.3. Mixing in D° -> K+tt tt° Decays 

The BABAR experiment has just reported further evidence for mixing in WS D° decays to K + ir~ir°. [7]. For 
this three-body channel, each decay is represented as a point in a Dalitz plot, defined by the squared invariant 
masses (sq,s + ) of the K + n~ and K + n° systems, respectively. The final state / and the value for the strong phase 
Sf = Sktttt + Arg (Af/Af) in Eq. (3) is then unique at each such point. The CF and DCS decay amplitudes Af and 
Af, respectively, can be determined from fits to their Dalitz plot densities, but the constant phase Sktt-k cannot. 

The BABAR collaboration identified a sample of - 3, 000 WS and - 660, 000 RS D° -> Ktttt events in a tight 
signal region with flavours identified by their origin from D* + decays. The RS sample was very clean, though the 
WS sample had a 50% background, mostly from "mis-tagged" RS decays. 

An isobar model for the amplitude Af was determined from a time-integrated fit to the distribution of events on 
the RS (mostly CF) Dalitz plot. A second, time-dependent fit to the WS Dalitz plot then determined parameters 
for the DCS amplitude Af and x" , y" (x and y rotated by the unknown phase Sktttt) summarized in Table I. 

Table I: Mixing parameters x" = xcosSkitit + ysinSKnn and y" — ycosSunn — x sin Saw from the BABAR analysis of 
D° -f K+-K--K decays. These are related to x and y by an undetermined strong phase rotation 5k n-n ■ 



2 BABAR measured the asymmetry in AY" = (y^p — y C p)/(ycp + Vcp) ~ 2r°A T /(r+ + r ). 



x 1 = x cos Sktt + y sin Sktt v' = V cos Sktt ~ x sin 5k-k 



2.2. Lifetimes for Decays to CP-even Final States 




(4) 
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Sample 


x" 


y" 


D° and D° 
D° only 
D° only 


[ 2.6ltg;g(stat.) ± 0.39(syst.)]% 
[ 2.53lo s3(stat.) ± 0.39(syst.)]% 
[ 3.55tg;£j(stat.) ± 0.65(syst.)]% 


[-0.06i° g(stat.) ± 0.34(syst.)]% 
[-0.05lo.s?(stat.) ± 0.50(syst.)]% 
[-0.54±°;*°(stat.) ± 0.41(syst.)]% 



The significance of the mixing signal in the fit to the combined D° and D° samples, ignoring CPV, was 3.2a. 
Separate fits to each sample show no evidence for CPV. 

3. OTHER MIXING MEASUREMENTS 

Two other channels have recently been analyzed in which the mixing signal was of less significance, but where 
mixing parameters were, nevertheless, measurable. 

3.1. D° -> K S 7r-7r + Decays 

Analysis of this channel, pioneered by the CLEO collaboration [22] is of great interest due to the presence of 
CP-eigenstates in the final state. These tie the phase of the amplitudes for D° and D° decay amplitudes together, 
so that Sktttv in Eq. (??) is zero. Thus, measurements of x, y and the two CPV parameters \q/p\ and <f>K s ir-w+ can 
be made free of unknown phases. 

CLEO's analysis used a e + e" sample and found only an upper limit for mixing parameters (—4.7 < x < 8.6)% 

and (—6.1 < y < 3.5)% at 95% CL. The Belle collaboration has repeated the analysis with a 540/6~ 1 sample (534K 
events) [4], observing mixing at the 2.4<r significance. The parameter values obtained are summarized in Table II. 
It is noteworthy that the value for y here is not consistent with the world average value for y CP = 1.072 ± 0.275% 
discussed in section 2.2. 

Table II: Mixing parameters obtained from the Belle analysis of D° — > Kstt~it + decays. The third error is that due to 
uncertainties in the model used to define the Dalitz plot amplitudes. 



X 




|q/p| 


Arg(q/p) 


(0.80±0.29±g;g?±g-.12)% 


(0.33±0.24l° ;f 2 l° ° 6 8 )% 


0.86l^ro'°o 6 3±0.08 


/ -..+16+5+2NO 
I ^-18-3-4/ 



There was no evidence for CPV since the magnitude of q/p was consistent with unity and its phase was zero. 

3.2. CP-odd Final States 

A new result from Belle [23]. has examined lifetime differences between the decays to the CP-odd final state 
D° — > 4>Ks with decays to the CP-even background and <j> sideband regions. A novel method was chosen over a 
time-dependent analysis of the whole KsK~K + Dalitz plot like that used for the Ksir~ir + channel which required 
a flavour-tagged sample of D°'s from D* decays. The sample size available was thereby enhanced almost threefold. 

Approximately 130K D° — > K~K + K$ decays were extracted from a 697/6 -1 sample of e + e~ collisions, with no 
attempt to distinguish D° from D°. Lifetimes t a and t b were measured, respectively, for events in the K~K + invari- 
ant mass regions A (containing the <j>) and B (narrow sidebands above and below the <j> where CP-even ao(980)Ks, 
fo(980)Ks and non-resonant decays were dominant). 

Decay amplitudes Ai,2 = A(s , s+) ± A(s , s-)]/y/2 were obtained from a fit [24, 25] to the distribution of events 
in the K~K + Ks Dalitz plot. Here, so, s_ and s+ were, respectively, K~K + , KsK~ and KsK + squared invariant 
mass combinations. Neglecting CPV, |^4i.2| 2 should have approximately exponential time-dependences with lifetimes 
t~i,2 = t/(1 ± yep), and their interference A\Ai should be zero when integrated over the s+ coordinate. 
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A value for y CP was extracted from the difference between (t a and t b ) and the "CP- ness", of each region, obtained 
by integration of these amplitudes. 

The result, y CP = (0.21 ± 0.63 ± 0.78 ± 0.01, where the third uncertainty accounts for uncertainties in the decay 
model used to describe the decay amplitude A, was included in Fig. 1(c). This result ranks third in precision for 
measurements of y CP . It is important that it agrees quite well with all other results from the CP-even modes. 

4. SEARCHES FOR CPV IN CHARM MESON DECAYS 

SM expectations for CP asymmetries in time-integrated decays in the charm sector are small. Singly-Cabibbo- 
Supprcsscd (SCS) decays arc of particular interest since these allow for asymmetries, at the ~ 0.1% level, from gluonic 
penguins [12, 26-28]. Observation of CPV at current sensitivities would be an indication for new physics. 

The asymmetry is defined as 

f = r(£> -> f) - T(D -> /) 

For D° asymmetries, there could be contributions from both direct CPV and from mixing. The latter is small since 
A T ~ 0.1%. 

Experimentally, measurement of asymmetries at the 10 -3 level are limited by uncertainties in asymmetries in 
the detection and reconstruction of particles of opposite charge. Also, for D° decays, efficiencies for D* + tagging 
cannot be assumed to be the same as that for D*~ . Forward-backward production asymmetries, resulting from Z°/7 
interference and higher order loops in the production of cc quarks, results in asymmetries in the distribution of D 
decay products in regions of varying efficiency in the detector. 

Calibration of these factors used to rely upon Monte Carlo simulated event (MC) studies, with questionable 
assumptions about charge-dependent interaction effects, resulting in systematic uncertainties in ^4 CP 's in the 1 — 5% 
range. In the past year, new insights in using data rather than MC have led to reduction of these uncertainties to 
the 0.2 - 0.4% range [29]. 

Large samples of CF D° — > K~tt + decays, selected from the BABAR 386/6- 1 data sample, were used to provide 
information on the efficiency asymmetries in small momentum and angular range. The underlying assumption was 
that the integrated production rates for D° and D° were the same. Forward-backward production asymmetries were 
dealt with by measuring event yields in forward and backward separately. 

Asymmetries obtained in this way, are summarized in Table III. Other, recent measurements are also included. 

Table III: Recent measurements of CP decay asymmetries xlCT 2 . The first uncertainty is statistical, the second is systematic. 
Data were used to estimate asymmetries in efficiency in all but ref. [30] . 



Mode 


BABAR x 10 2 


Belle xl0~ 2 


CLEO xl0~ 2 


D° -> K K+ 

D°^7r-7r+ 
D°^K K+7T 
D^tt-tt+tt 
D+ -> K K+7T+ 


0.00 ±0.34 ±0.13 [29] 
-0.24 ± 0.52 ± 0.22 [29] 

1.00 ± 1.67 ±0.25 [32] 
-0.31 ±0.41 ±0.17 [32] 


-0.43 ±0.30 ±0.11 [31] 
0.43 ±0.52 ±0.12 [31] 

-0.43 ±0.41 ± 1.23 [30] 


-0.03 ± 0.84 ± 0.29 [33] 



These asymmetries are consistent with zero, with most systematic uncertainties much less than 1%. Of particular 
interest is the observation that both statistical and systematic uncertainties should scale with the square root of 
the number of events. Future measurements from flavour factories may well detect asymmetries at or above the SM 
limit. 
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5. SUMMARY 

The Heavy Flavour Averaging Group (HFAG) has combined a wide range of mixing observables (28 in all) including 
all those discussed here. To do this, they made a x 2 fit [34] to obtain values for the underlying mixing parameters 
(x, y, \q/p\,etc.) that best described these. The results are best summarized in the \ 2 countour plots in Fig. 2 for 
the fit that allowed for CPV . The combined evidence for mixing is now compelling. The probability for x = y = 




Figure 2: \ 2 contours for fit to mixing observables. (a) x vs. y and (b) |g/p|vs. Arg(q/p). This figure is taken from the Heavy 
Flavour Averaging Group (HFAG) [19] 

(no mixing) is excluded from a fit with no CPVdX the 9.8a level. The present data are also consistent with the 
absence of CPV. No evidence for CPV in time-integrated decay rates for neutral D mesons exists either, though the 
measurements are approaching the interesting level of sensitivity that could confront SM estimates. Furthermore, 
systematic as well as statistical uncertainties are expected to shrink as more data at present and future flavour 
factories is accrued. 

We should take this as a need to continue to press on with these studies in case there is yet a hole in the SM at 
the energy scale that can be probed with results from a flavour factory. 
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